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We provide convincing empirical evidence that long range interactions strongly enhance the rec¬ 
tification effect which takes place in mass graded systems. Even more importantly the rectification 
does not decrease with the increase of the system size. Large rectification is obtained also for the 
equal mass case and with graded on-site potential. These results allow to overcome current limi¬ 
tations of the rectification mechanism and open the way for a realistic implementation of efficient 
thermal diodes. 
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The great development of modern electronics with the 
invention of electric diodes, transistors and other nonlin¬ 
ear solid-state devices has led to impressive scientific and 
technological achievements, with strong impact in our 
daily lives. However its counterpart, aimed at the ma¬ 
nipulation of the heat current and phononics, more than 
a decade after the first proposal of the thermal diode flj, 
remains at a near-standstill in terms of management of 
thermal circuits, basically due to the lacking of efficient 
and feasible diodes. 

The thermal diode is the basic tool for the manipu¬ 
lation and control of the heat flow. In this device heat 
flows preferably in one direction, i.e., the magnitude of 
the heat current changes as we invert the sample between 
two thermal baths. The most recurrent proposals of ther¬ 
mal diodes are based on the sequential coupling of two 
or three segments with anharmonic interactions jlH3( and 
they have to face several problems which are difficult to 
overcome. In particular their rectification power is small 
and rapidly decays to zero as the system size increases. 
For this reason, more and more efforts have been devoted 
to the search of alternative systems or rectifying mecha¬ 
nisms, however without remarkable results. 

In this letter, we suggest a possible way to overcome 
this conundrum, namely we identify general and simple 
ingredients to build a diode with great rectification effect 
which moreover does not decay to zero with increasing 
the system size. By considering one dimensional anhar¬ 
monic chains of oscillators, we provide empirical evidence 
that graded mass distribution and long range interparti¬ 
cle interactions, lead to a substantial improvement of the 
thermal rectification phenomenon. It is interesting that a 
comparable very large rectification factor is also obtained 
if asymmetry in the system is introduced in a different 
way, e.g., by considering the equal mass case with graded 
on-site potential. 

We would like to stress the following points. First, 
since the pioneering work of Debye, anharmonic chains 


of oscillators are the most basic and recurrent models 
for the study of heat conduction in insulating crystals as 
they are believed to mimic, to some extent, the behavior 
of real materials. 


Second, we recall that graded materials, i.e., inhomo¬ 
geneous systems whose structure changes gradually in 
space, are not only abundant in nature, but they can 
also be manufactured. They have attracted considerable 
interest in many areas: material science, engineering, op¬ 
tics, etc. 0. In particular, a simple type of graded diode 
has been already experimentally built Q: it consists of 
a carbon and boron nitride nanotube, inhomogeneously 
mass-loaded with heavy molecules. Unfortunately, the 
rectification observed there is very small. 


Finally, we emphasize that long range interactions are 
not uncommon in physics: they are at the origin of inter¬ 
esting properties and phenomena in classical and quan¬ 
tum physics, in equilibrium and nonequilibrium phase 
transitions, etc. They are important at all scales: from 
cosmology [6] to nanoscience [3 and quantum informa¬ 
tion 0. 0 . More specifically, in material science, present 
technology allows to fabricate materials with different 
properties and even different types of long range inter¬ 
actions. A noticeable example is given by some Coulomb 
crystals which are systems intensively investigated and 
also proposed as candidates for a many-qubit quantum 
processor. In particular the engineered Coulomb crystal 
of ref. [lot , given by a trapped beryllium ion system, 
exhibits Ising like interactions with ,/,j oc 1 /rfj, where 
rij is the distance between spin pairs and 0 < a < 3. 
Long-range magnetic dipolar interactions are present in 
the Ising pyrochlore magnets Dy2Ti20i and H 02 T 12 O- 1 , 
which are also of current interest due to the so called 
spin ice behavior 11]. An additional example is provided 
by nanomagnets such as Permalloy, materials which can 
be lithographically manipulated jl2| in order to display 
frustrated arrays and exotic states. In such systems the 
interaction typically decays polynomially as 1 /r 3 , where 
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FIG. 1 : (Color online) The schematic plot of our model. At 
the left and right ends of the chain there are two leads with 
particles of the same mass mi (on-site potential strength 71) 
and mjv (71V), respectively. Particles in the leads have only 
nearest neighbor interaction. We have verified that in actual 
computations it is sufficient to have one particle in each lead. 


r is the distance between the nanodisks. 

The system we will discuss here is a one-dimensional 
chain of N oscillators described by the Hamiltonian 


(gj - g i) 

2 + 2\i — j\ x 


(1) 


where qj is the displacement from the equilibrium posi¬ 
tion of the jtli particle with mass rrij and momentum 
Pj. The exponent A gives the decay power of the inter¬ 
particle interaction, while 7 j measures the strength of 
the on-site potential [l3(. The system of anharmonic os¬ 
cillators is connected to two leads, whose particles have 
nearest neighbor interactions only so that the heat flux 
through the leads has the usual simple definition fl^-fl 6 i| . 
Two Langevin heat baths at different temperatures 
Tl and Th are attached to the two ends of the leads; we 
take T l = T{ 1 - and T H = T( 1 + (See Fig.l). 

The Hamiltonian of the lead reads Hi ea d = 
S"=i (27IF + + k( ' qi+1 ~ qi ' > ^, where n is the number 

of particles in the lead and k is the nearest neighbor 
coupling strength. The main system is coupled to each 
lead by a harmonic spring with the same nearest neigh¬ 
bor coupling strength k. In the simulations, k = 1 is 
adopted. We have checked that the main results of this 
paper do not depend on the value of k and on the number 
of particles in the lead. 

In view of previous results [2] , it is expected that in a 
system with graded mass distribution, e.g. m 1 < m 2 < 
... < to, y, thermal rectification will be present, even 
for the simple case of nearest neighbor interaction (NN). 
Now, with long range interactions (LRI), we conjecture 
that the presence of new links (interactions) among differ¬ 
ent sites creates new channels for the heat transport, that 
favor the heat flow. Besides that, in a graded system, the 
new channels connect distant particles with very different 
masses. Therefore new, asymmetric channels, are created 
which in turn favors the asymmetric flow, i.e., rectifica¬ 
tion. Hence, by introducing long range interactions in a 
graded system, we expect an increase of the thermal rec¬ 
tification. Moreover, as we increase the system size, new 
particles are introduced that, in the case of long range in¬ 
teractions, create new channels for the heat current which 


may avoid the usual decay of rectification with increas¬ 
ing the system length. In conclusion, systems with long 
range interactions may represent a relevant step on the 
way of constructing efficient thermal diodes. 

We would like to remark that the possibility of rec¬ 
tification growth with long rang e interactions has been 
suggested in a recent work [19( where, due to analyti¬ 
cal difficulties, only an asymptotic study of some simple 
model has been presented (chains of oscillators with self- 
consistent inner reservoirs). Moreover, only infinitesimal 
differences of temperature are considered, leading to in¬ 
finitesimal rectification. Therefore the role of long range 
interaction in the rectification phenomenon remains to 
be understood. 

Let us now turn to numerical simulations of our model 
in order to determine its thermal transport properties 
as a function of different parameters such as temper¬ 
ature gradient, mass gradient, etc. In particular we 
will compute the rectification factor, defined as f r = 
- J+ J J x 100%, where J + and J_ represent, respec¬ 
tively, the larger heat flow and the smaller heat flow, 
which are measured by inverting the temperatures of the 
baths at the two ends of the chain. We found that the 
larger heat flow J + is obtained when the heat bath at 
higher temperature is on the side of the heavier mass. 

In our simulations the system is evolved for a long 
enough time to ensure that it has reached the stationary 
state. After that the heat current is measured as the av¬ 
erage heat flow through the two leads. More specifically, 
we verified that the average heat flow into the system 
is equal to that out of the system. We have also per¬ 
formed a series of numerical test in order to ensure that 
our results do not depend on numerical technical details 
or on the choice of the heat baths (for example we have 
compared with simulations using Nose-Hoover heat baths 
and obtained consistent results). In our simulations, free 
boundary conditions and velocity-Verlet algorithm are 
used. We have compared with simulations using fixed 
boundary conditions, and with simulations using Runge- 
Kutta algorithm of seventh to eighth order, and obtained 
consistent results. 

As commonly adopted in the literature, we use here 
the dimensionless units. 

We start our numerical studies by investigating the role 
played by the decay power A of the interparticle interac¬ 
tion. We expect that, as the decay power is increased 
and, correspondingly, the interaction between sites be¬ 
comes smaller, the rectification factor should decay. Such 
scenario is confirmed in Fig.2 where we plot the rectifica¬ 
tion factor as a function on the decay exponent A. Notice 
however that for A close to one, the rectification factor 
reaches a maximum value and decreases again when A 
approaches one. The reason of such behaviour is not 
clear to us. We surmise that this behaviour is connected 
to the finite system size and indeed the position of the 














3 



A 


FIG. 2: (Color online) Rectification factor versus the decay 
exponent A of the interparticle potential. Here T = 5, = 

0.9, mi = 1 and 7 j = 1 for all j. Squares are for N = 64, 
mjv = 10; triangles are for N = 256, miv = 5; stars are for 
N = 128, m jv = 3; and circles are for N = 64, mjv = 2. 
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FIG. 3: (Color online) Dependence of the rectification factor 
on the temperature gradient. Here T = 5, Tl = T(1 — 

Th = T( 1 + ^f-), N = 64, mi = 1, mjv = 10. 7 j = 1 for all 
j. Triangles are for LRI with A = 1.2, circles are for LRI with 
A = 1.4, squares are for the NN case. 

maximum becomes closer to one as we increase the sys¬ 
tem size. 

We now analyze how the thermal properties depend on 
the main parameters which determines the asymmetry in 
the system. In particular we will compare the behavior 
of the thermal rectification factor for the LRI-chain and 
the corresponding chain with nearest neighbor interac¬ 
tion only. The dependence of f r versus the temperature 
difference is depicted in Fig.3 for a system with size and 
mass gradient fixed. The f r amplification due to the 
presence of LRI is remarkable. 

As argued before, as we increase the asymmetry in 
the system, we expect to increase the rectification. The 



FIG. 4: (Color online) Dependence of the rectification factor 
on the mass gradient. Here mjv = 10 , N = 64, T = 5, 
= 0.9 and 7 j = 1 for all j. Triangles are for LRI with 
A = 1.2, circles are for LRI with A = 1.4, squares are for the 
NN case. 

effect is clearly seen in Fig.4, where the mass gradient is 
increased by decreasing mi with fixed tun- 

All the above results make transparent the significant 
gain in the thermal rectification due to long range inter¬ 
actions: for example for the case Th = 9.9, Tl = 0.1, 
N = 64, mi = 1, mu = 10, A = 1.2, we get the quite 
large value f r = 4638% f2c|. 

We now turn to the investigation of rectification prop¬ 
erties as a function of the system size. This is particularly 
relevant since an unpleasant feature of previous models 
of rectifiers is the decay of rectification with the increase 
of the chain length. The results are shown in Fig.5 in 
which the mass gradient is fixed. We see that the pres¬ 
ence of LRI prevents the decay of f r with the system 
size. Strictly speaking we cannot make any claim for 
larger system sizes. However it is clear from Fig.5 that 
the N-dependence for the LRI case is qualitatively differ¬ 
ent from the NN case where the decay of the rectification 
factor with N is immediately evident. 

We would like to stress that the rectification effect de¬ 
scribed here is caused by asymmetry and therefore it 
should be present in different asymmetric chains. For ex¬ 
ample, the rectification phenomenon should be present 
also in systems in which the asymmetry in the parti¬ 
cle mass distribution is replaced by an asymmetry in 
the on-site anliarmonic potential. In order to check 
this prediction we have considered our model ( 1 ) with 
equal masses and with graded on-site potential, namely 
we took 7 i linearly increasing with i. More precisely 
7i = 1 < 7i < 7 max =7 N(i = 1,2,..., N). The results 
are shown in Fig .6 which shows a strong rectification ef¬ 
fect. This result is particularly interesting since it shows 
that rectification may be implemented in a symmetric 
chain in which asymmetry is induced by a transverse ex- 
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FIG. 5: (Color online) Dependence of rectification factor on 
the system size N. Here T = 5, = 0.9, mi = 1 and 7 j = 1 

for all j. Triangles are for LRI with A = 1.2, circles are for 
LRI with A = 1.4, squares are for the NN case. The mass 
gradient is 1/64. (mjv = 2 for N=64; rriN = 3 for N=128; 
rriN = 4 for N=192; rriN = 5 for N=256.). 
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FIG. 6 : (Color online) Rectification factor versus on-site po¬ 
tential coefficient 'ym.ax, for equal-mass system. Size N = 64, 
T = 5, = 0.9, 7 min = 1. (From particle 1 to particle N, 

7 increases linearly. When the hot heat bath is on the 'Ymin 
side, the heat flux is larger.) 


ternal field which is graded along the chain. 

In this paper we have considered a one dimensional 
model for reasons of numerical convenience. However 
the rectification effect should be present in larger dimen¬ 
sional systems since asymmetry, which is at the root of 
the rectification mechanisms, can be realized also in two 
or three dimensions. In this connection we stress that our 
results are obtained for graded chains of anharmonic os¬ 
cillators without any particular requirement or technical 
restriction. 

The main impasse in the development of a more intense 
experimental research activity on heat current manipu¬ 


lation is the low rectification factor which characterizes 
current proposals of thermal diodes as well as the fact 
that the rectification factor decreases with increasing the 
system size. In this paper we propose a way to overcome 
this difficulty. More precisely, we show that the presence 
of long range interactions in a graded mass (or graded 
structured) system can considerably increase the recti¬ 
fication factor and avoid its decay as the chain length 
grows. Finally, we stress that the proposed mechanism is 
not restricted by geometrical considerations or by space 
dimensionality. We are confident that our results will 
stimulate experimental research on this subject. 
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